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I»  Introduction 

In  an  earlier  paper  [2]  we  showed  how  the  notion  of 
Taylor  Instability  can  be  used  to  explain  the  break-up  of 
accelerated  thin  liquid  sheets  into  drops,  and  how  on  the 
basis  of  this  theory  the  drop  sizes  can  be  estimated.   The 
idea  underlying  the  calculation  is  as  follows:   One  considers 
a  plane  layer  accelerated  in  a  direction  normal  to  its 
surface  by  imposing,  e«g.,  a  pressure  difference  on  opposite 
surfaces.   The  zero  order  motion  —  which  incidentally  Is 
an  exact  solution  for  the  plane  layer  —  Is  a  parallel  flow 
with  the  velocity  of  bounding  surfaces.   It  is  next  argued 
that  the  flow  actually  deviates  from  the  zero  order  solution 
because  either  the  bounding  surfaces  are  not  perfectly  plane, 
or  the  pressure  on  the  boundary  is  not  exactly  constant,  or 
because  of  some  random  perturbation  that  may  occur  at  the 
outset  or  during  the  motion.   To  see  what  happens  to  the 
bounding  surfaces,  one  considers  the  first  order  perturbation 
which  satisfies  linear  equations  and  is  represented  by  a 
series  (or  integral)  of  normal  modes.   Some  of  these  modes 
are  found  to  be  unstable  in  the  sense  th-nt  their  amplitudes 
grow  unrestrictedly  with  time.   Among  tries©  one  or  more  may 
be  called  most  unstable  in  the  sense  that  they  grow  most 
rapidly,  VThen   the  acceleration  of  the  layer  is  constant  the 
unstable  modes  grow  exponentially  and  the  most  unstable  modes 
have  the  largest  exponent.   The  most  unstable  modes  are 
assumed  responsible  for  break-up  of  the  sheet  by  pinching 
it  into  pieces  which  eventually  become  spherical  drops  under 
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2. 

the  Influence  of  surface  tension. 

On  the  basis  of  the  above  mechanism  one   gets  an  estimate 
of  the  number  of  drops  produced  per  unit  area  and  subsequently 
a  formula  for  their  radii.   This  formula  is 

/    ?    1/3 

where 

h  —  thickness  of  the  sheet 

1>  —  coefficient  of  surface  tension 
(Pi-Pp)  —  difference  of   pressures  exerted 
on  different  sides  of  the  layer. 
The  formula  (1)  was  obtained  under  the  assumption  th'jt  the 
liquid  is  incompressible  and  invlscld.   As  It  indicates,  the 
notion  of  surface  tension  plays  a  basic  role,   V/lthout  It, 
Indeed,  the  notion  of  maximum  instability  vjould  not  exist. 
Concerning  this  see  also  [1]. 

A  formula  such  os  (  )  is  of  considerable  practical 
Importance,  yet  its  validity  is  in  doubt.   Even  If  one  does 
not  question  the  correctness  of  the  mechanism  used  to  explain 
the  process  of  dropletizatlon,  a  doubt  arises  from  an  un- 
qualified extrapolation  of  the  notion  of  instability  by  means 
of  which  the  formula  has  been  obtained.   Although  the  exis- 
tence of  instabilities  is  unquestionable,  it  is  not  at  all 
certain  that  small  disturbances  will  grow  sufficiently  to 
produce  pinching  as  the  linearized  perturbation  theory 
Indicates.   That  an  indiscriminate  growth  of  perturbations 
may  be  arrested  is  indicated  by  the  experiments  of  Levis 
[5].   Since  then,  Layzer  [I4.]  has  shovm  by  means  of  a  more 
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refined  theoretical  considerations  that  the  growth  of  iin- 
stable  modes  is  Indeed  arrested.   The  calculations  by 
Pennington  and  others  [6]  point  to  the  same.   Layzer's  and 
Pennington's  papers  deal  with  the  Taylor  case  —  half 
infinite  liquid  medium  without  consideration  of  effects  of 
surface  tensions  —  and  their  result  Is  thus  much  more 
striking. 

Nevertheless,  accelerated  sheets  do  desintegrate. 
The  outstanding  question  is  then  whether  the  mechanism 
of  desintegratlon  proposed  In  [2]  Is  justifiable;  and  whether 
the  formula  (1)  has  at  least  a  limited  practical  value.   In 
[2]   a  precaution  is  taken  by  stating  that  the  treatment 
applies  to  thin  sheets.   What  is  needed  is  a  qualification 
of  what  this  means  in  this  problem  and  a  description  of 
what  happens  when  the  sheet  is  not  thin.   In  this  paper  we 
shall  present  a  partial  answer  to  these  questions. 

While  formula  (1)  predicts  a  natural  desintegratlon 
of  sheets  into  drops  of  definite  size,  a  question  arises 
about  what  happens  to  an  accelerated  sheet  when  initially 
one  surface  of  the  layer  is  deformed  in  a  controlled  way, 
by  introducing,  e.g.,  slight  dimples  distributed  periodically. 
This  poses  an  initial  value  free  boundary  problem,  which  is, 
of  course,  very  difficult  to  handle.   VJe  consider  here  a 
very  special  problem  of  this  type,  namely  that  In  which  one 
surface  is  kept  rigid,  while  the  other  surface  has  a 
sinusoidal  profile.   That  is,  the  profile  is  given  initially 
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by   the   equation 

(2)  2  =   a  cos  -^      . 

We  can   treat   also  the   case   vhenthe   equation  of  the 
surface   Is   Initially 

(3)  z  =    a  cos  —  cos  -^   , 

but  this  vould  lead  to  technical  complications  without 
adding  to  the  understanding  of  the  behavior  of  the  layer. 

It  turns  out  that  the  equation  for  the  surface  at 
any  time  can  be  expressed  as  a  f  omial  power  series  in  a/1 
vhlch  probably  has  an  asymptotic  charactero   Furthermore, 
with  the  assumed  Initial  data  (2)  (or  (3),  vhich  we  do  not 
consider  here),  and  in  these  cases  only,  this  series  can  be 
computed  to  any  order.   We  carried  our  c considerations  up 
to  terms  of  third  order,  and  to  this  order  our  results  seem 
to  a gree  with  o ther  theoretical  considerations  and  with 
some  experimental  results,   (It  is  not  known  ih  ether  the 
consideration  of  higher  approximations  would  improve  the 
accuracy  of  these  results.) 

In  Section  2  the  mathematical  problem  is  formulated. 
In  Section  3  we  consider  the  instability  theory,  i.e. 
the  first  order  approximation  as  applied  to  the  case  of 
surfaces  having  one  dimensional  profile.   In  Section  l^   we 
present  the  method  of  solution  t^Jilch  is  followed  in  Section 

5  with  explicit  results  up  to  terms  of  third  order.   Section 

6  contains  a  discussion  of  these  results. 
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We  consider  sheets  both  of  finite  and  Infinite  thickness. 
As  to  the  latter,  our  results  agree  satisfactorily  with  those 
of  Pennington  [5]  who  used  a  high  speed  electronic  computer 
to  study  the  growth  of  instability.   They  also  show  an  arrest 
of  instabilities  as  predicted  In  [5]  and  [6]. 

Concerning  sheets  of  finite  thickness,  our  study  of 
their  time  history  shows  at  first  an  amplification  of  their 
initial  deformation  followed  in  many  cases  by  the  formation 
of  Jets  at  the  trough  of  the  layer,  see  Figures  8,  10,  and  12, 
This  feature  is  very  plausible,  as  can  be  seen  from  the 
following  qualitative  description  of  the  layer  based  on  t he 
assumption  of  incomprensibillty  of  the  liquid,  see  Fig.  1-a, 
«-b,-o,  page  6, 

Figure  la  shows  the  initial  configuration.   The  liquid 
Is  at  rest  but  the  lower  surface  is  deformed.   The  upper 
surface  will  be  kept  plane  during  the  motion.  As  the  layer 
is  accelerated  upwards,  the  deformations  of  the  lower  (free) 
surface  will  grow,  due  to  instability.    At  least  during  an 
early  stage  of  the  motion,  liquid  will  flow  generally  away 
from  the  troughs  into  the  crests  of  the  surface*  This 
Ls  indicated  in  Plg^  lb  by  arrows  along  the  streamlines.  As 
the  motion  progresses,  either  the  sheet  will  break  or  the 
growth  of  the  instabilities  will  be  arrested.   The  first 
possibility  is  exactly  the  basis  for  the  theory  of  break-up 
proposed  in  [2].   In  the  other  case  there  will  be  a  time  at 
which  the  velocities  in  the  fluid  vl  11  reverse  their  direction. 
At  that  time  due  to  high  velocities  at  the  troughs  of  the 
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7. 

layer,  jets  will  be  produced  as  Indicated  in  Pig.  Ic .   A 
similar  situation  was  observed  experimentally  by  T.  Holland^  ' 
vho  used  high  speed  photography  to  study  the  early  history  of 
accelerated  sheets. 

As  to  the  later  history  of  the  sheet,  our  results  seem 
to  confirm  to  a  certain  degree  those  of  the  instability 
theory.  As  time  progresses  the  deformation  becomes  so  large 
that  the  sheet  breaks  into  filaments  by  pinching.   If  one 
starts  with  a  deformation  xhose  wavelength  Is  approximately 
equal  to  the  critical  wavelength  given  by  the  instability 
theory,  the  breakdoi-m  is  into  approximately  one  filament  per 
wavelength.   If  it  is  twice  the  critical  wavelength,  the 
breakdown  is  into  two  filaments,  which  are  not,  however,  of 
same  size.   If  finally,  it  is  three  times  the  critical  wave- 
length, the  breakdown  is  into  three  filaments.   This  is  In 
particular  true  when  the  sheet  thickness  Is  less  than  the 
critical  wave  length.   For  thick  sheets  our  theory  may  not 
be  adequate.   In  all  cases  our  result  shows  a  slower  rate  of 
deformation  than  that  predicted  by  the  instability  theory. 

Although  we  have  treated  only  the  case  of  a  one  dimen- 
sional profile,  which  leads  to  a  break  of  sheets  into  fila- 
ments, our  results  indicate  that  they  apply  equally  well  to 
the  two  dimensional  profile,  when  a  deslntegratlon  into  drops 
is  expected.   Our  conclusion  is  thus  that  the  formula  (1)  is 
adequate  for  the  case  when  the  sheet  thickness  Is  small 
compared  to  the  critical  wavelength. 


(1)  Oral  commtmloatlon.   The  author  is  indebted  to  Dr.  T. 
Holland  for  making  these  results  available. 
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2.  Formulation  of  the  free  boundary  problem, 

zA  Consider  a  layer  of  incom- 

pressible, inviscid  liquid  contained 

hi  between  the  surfaces  P,  and  P„t 

p  ±  d 


Liquid         2 


:r  ■  -  •  -^^-x r— X  -~>  P, 


Assume  that  the  layer  is  accelerated 

I  upwards  with  the  constant  accelera- 

i  _       tion  a  >  0.   The  upper  boiJindary  is 

1   kept  plane,  while  on  the  lower,  free 

boundary,  there  is  exerted  a  fixed 
dirdction 
of  arrow  pressure  p  augmented  by  pressure 

due  to  surface  tension.   The 

^^S'    2  equations  of  the  boundaries  are 

1  2 

(2.1)  z  s  h  +  -^^t     upper  boundary  Pp 

1      2 

(2.2)  2  =  -^at     +   Z(x,y,t)    lower  boundary  P,    , 

where   h   is    the    average   thickness  of  the    layer. 

The    object   is    to   find   Z(x,y,t)    for   all   t  given  Z,Z.    at   t  =   0    . 

In   the   liquid,    the    velocity  field   is   derived   from  a 
potential    cj)  satisfying  Laplace's    equation, 

(2.3)  A  ^  =    0« 

The  kinematic  boundary  conditions  are 

(2.[|.)  4^  =  at    on  F^    , 

(2.5)  i^-Z^i^^Zycty  =  at+Z^   on  P^  . 

On  the  free  boundary  we  have  furthermore  the  dynamic  condition 
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stating  continuity   of  pressure, 

(2.6)  f>{^^+^{V^)^  +  p   +  2TH  =   f(t),      on  P^    . 

Here  T   la   the    coefficient   of  surface  tension  and  f(t) 

is    an    arbitrary  function  of  time   whose   choice   Is    determined 

by  the  normalization   of   the   potential    (j),      H  is    the  mean 

curvature   of  the    free  surface   given  by 

.3 

(2.7)  H  =  |((l^z2)Z^^.2Z^ZyZ^^y-H(l^z2)Z^)(l.z2.z2)    ' 

We   restrict   the   study   to   the    special  case    when 

(2.8)  Z(x,y,o)    =  X  cos   ^   ,    Z^(x,y,o)    =   0   . 

Since  y  does   not   appear  explicitly   In   any   of  the   above 

equations,   the  solution  will  be   Independent  of  y.      H  will. 

In  particular,   take    on  the   simpler  form 

.3 

(2.7')  H=|(l.z2)      Z^^    . 

It   Is   convenient   to   consider  the  problem  in  a  frame 
moving  with   the    upper   surface,    and    to    introduce  dimenslonless 
varia>>les.      Chose  5/    and  T  =  /o/a      for   units   of  length  and   time 
respectively.      Let 

x  =   H,  y  =  9y,    z  =  iz  +  iat^t^,   t  =     t 

,2   -  _ 

(2.9)  ci  =  -^  (l>(x,y,z,t)    +  ^ratz  ^•P(t) 

Z  =    JvZ(x,y,t). 

The  conditions  (k)  tiS) ,    and  (6)  become,  on  dropping  all  bars 
and  setting  for  convenience 

F'(t)  =|a2',2;2  _  £  ^  l^(^^j  ^ 
2         |-    f 
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i2.sn 

(2.6M 

where 

(2.10) 


^     =  0,   for   z  =   V      , 


^ 


Z  (J)      -  z  ^     =   Z.     ,   for   z  =    Z      , 


z  x"x         y^y  t 


4^  +  i   (  X^  <J,)2  +   2  +  2kH  =   0,   for   2  =   Z   , 


^=1 


< 


k  = 


L 


pa 


F 


Equations    (2,3)    and    (2.7}    (2.7')    remain  unchanged. 
The   special    initial  conditions    (2,8)    become 


a 


(2.8n      Z(x^y,0)  =    e  cos  x,   Z^(x,y,0)  =    0  ivihere   e  =  ^ 

The  problem  contains   three   parameters,      v   Is    the   thick- 
ness  parameter  measured   in  units   of  the    vave    length   of   the 
initial  disturbance,   e   is   the   amplitude    of  d isturbance,    and 
k  may  be   conveniently  called  the   instability  parameter. 
3»   Instability  theory 

If  ve   disregard   the  initial   conditions    (2e8'),    (|)  =  0, 
Z  3  0  give   a  possible   motion  o f  the    layer  -  the  so-called 
zero  order  motion.     The  instability  theory  studies    the 
possible   deviations   from  the   zero  order  solution  under   the 
assinnptlon  that  these   deviations   are    so  small  that  all 
equations  may  be    justifiably  linearized*     Assume   that 
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4  =  ;      g(t)    cosh  X(v-z)\j/.  (x^y)      , 


z  =  II^x^*^    ^x^^'^^      ' 


^ere  X  Is  some  Indexing  system.  Equation  {2^^})   is  satisfied 
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and  (2,3)  will  be  satisfied  if  ^     Is  a  solution  of 

(3.2)  \l/^^+       I'yy       +       X^\l/       =       C        . 

The  remaining  tvjo  boiindary  conditions  now  imply,  on  linear- 
ization that 

(3.3)  -Xg  3inh  XV  =  f  , 

(3.U)  g  cosh  XV  +  (l-kX^)f  =  0 

Solving  for  f,  one  gets 

f(t)  =  fol«°**^o2«'°* 

where        p  P 

a  =  X  tanh  Xv(l-kX  )  . 

2   -1 
VFhen  X  <k    ,  a  will  be  real  and  therefore  the  time 

factors  corresponding  to  those  values  of  X  will  grow  exponen- 
tially. The  most  unstable  mode  will  correspond  to  that  value 

of  X  which  maximizes  a,  i.e.  for  which  -r—(<^    )='0'.   X   is 
m  dX  rn 

the  solution  of 

(3.6)      (l-3kX^)  th  XV  +  vX(l-kX^)(l-th^Xv)  =  0, 

When  X   V   is    large,    we  have    approximately 
Set 


Substituting    In    (6),    one   gets,    after  some   transformations 

2       cr  +  sinh  a 


(3.9)  u^  =  3- 


rrU  +   sinh  C7 


"  '■ 


i;    i    .r;U. 


t  '"/■  <- 1. ;    *• 
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where 


a  = 


2v^ 


Clearly  then  1  </^<    /-^   ,  so  that 
(3.10) 


—  <  ^  <  — 

For  V  =  CD  (Taylor  case),/-t.=  i,   Llkevlae,  for  k  =  0  (i.e. 
vhen  T  =  0),/^=  1.  (In  the  latter  case  the  notion  of  maxi- 
mum instability  does  not  exist.)   The  table  belov;  gives  the 
correction  factors/^  for  those  values  of  v  and  k  which  are 
used  for  illustrative  purposes  in  further  discussion, 
Section  6. 


Table 
m     J  3k 

(3.1) 

/ 

I 

k 

V  =    1 

V  =    .2 

0 

CO 

1 

1 

1/27 

3 

1.003 

1.177 

1/12 

2 

1.0U5 

1.203 

1/3 

1 

1,12 

1.219 

2/3 

.71 

1.16 

1.222 

1 

.58 

1.18 

1.223 

00 

0 

/I 

n 

V  2 

Assume  now  that  the  lower  surface  is  a  cylinder  with 
generators  parallel  to  y  axis,  The  profile  of  the  surface  is 
then  one  dimensional,  and  the  liquid  flow  -  two  dimensional. 
Prom  the  equation  (2)  we  get,  assuming  the  normalization 
\Ko)  =  1,  \Jr^o)  =  0  , 
(3.11)  ^x  ~  °°^  ^^ 

According  to  the  instability  theory,  the  sheet  id  11  break 

p-ir 

into  filaments, of  width  ^  ,  that  is,  into  X  filaments  per 

X^    '  '  m  ^ 

m 

section  of  width  2ii, 
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Ij..-  Method  of  'Solution. 

We  now  turn  back  to  the  free  boundary  problem  formulated 
in  Section  2.  Assume  that  Z  and  (J)  can  be  expanded  In  formal 
power  series  (convergent  or  asymptotic)  in  e  and  that  further- 
more the  coefficients  in  these  expansions  are  themselves 
expanded  in  Fourier  series  in  x,  vlth  period  2  it  •  Thus,  quite 

generally, 

05     oo 
(l4.,l)      Z(x,t)  =  ^  e  T[  (a   (t)  cos  mx  +  a^^(t)  sin  mx) 

n=0   m=0 

CD      CD 

(1+.2)   4(x,z,t)  =  H  e"T"  (b„^(t)  cos  mx 

m='0   m=C 

+b (t)  eln  mx)  cosh  m(v-e), 

nm  '  ' 

Each  term  in  (i;.2)  satisfies  already  the  oondltlons  (2»2,2tl\.^  ) , 
The  problem  thus  is  reduced  to  the  determination  of  the  a's 
»  and  b's  in  such  a  way  that  the  remaining  boundary  conditions 
(2,5',  2,6')  be  also  satisfied.   The  initial  conditions  (2,8') 
Imply  that  at  t  =  0 : 

(l4.,3)   a   =a   =:a   =a   =0,  for  all  n  and  m  except  that 
nm    nm    nm    nm    ' 

To  determine  the  a's  and  b's,  vje  substitute  (i^-.l)  and 
(I|.,2)  in  (2,5',  2,6'),G:arry  out  the  intended  operations, 
group  terms  of  like  powers  of  e,  and  lastly  expand  each 
coefficient  in  a  Fourier  series  in  x.   This  results  in  each 
case  in  an  equation  nf  the  form 

00     oo 

n=o   m=« 


:  -IT};; 


a  .-• 


Hno     i-? 


u     .o:     i--f 


z  -r. 


•J     .  ^    f .  - . 


'.  y.': 


Zb::-^i 


.■J.r 


'^  .  *-  ••*       •  . 


Ill     ;jf; 


I  • 


•     -\ 


.*C.      > 


.1  '.'i 


'        ■■•  ' 


•  •:»•  * 


T.There  the  X's  depend  on  the  a'a  and  b's  only.   Equation 
(I4..5)  Implies  that  all  the  X's  are  Identically  zero,  thus 
leading  to  an  Infinite  set  of  equations  for  the  a's  and  b's. 

On  setting  6=0,  one  gets  the  case  when  t he  Initial 
free  surface  Is  undeformed.  In  that  case  the  solution  la 
the  zero  order  solution, 

Z  =  0,  i  =  0   ; 

hence 

(li,6)        a   sa   =h   =b   =0. 
^^   '         on  -  om    om    om 

With  an  arbitrary  Initial  Z  there  would  be  for  every  n  >  1 
an  Infinite  number  of  equations  to  solve.   Thus  In  general 
we  would  not  be  able  to  obtain  any  results  beyond  the  terms 
of  first  order  in  e.   The  choice  of  the  particular  initial  Z 
given  by  equation  (2,8')  Is  motivated  by  the  fact  that  in 
this  case  (and  in  this  case  only)  for  each  value  of  the  first 
index  there  is  only  a  finite  number  of  not  identically 
vanishing  a's  and  b'sj  furthermore,  the  equations  for  these 
can  be  solved  in  succession.   More  precisely  the  following 

is  true: 

l,a       =b       =0,   for  all  in  and  n. 
nm  nm  ' 

2,  a        =0,    for  all  n 

no  ' 

3.  The  pairs    (a      ,b     )  satisfy  equations    of   the 
^  ^  ^    nm*    run  *'      ^ 

form 


(k.7) 


a        +    (m   sinh  m   v)b       =   f  _, 
nm        ^  '    nm         nm' 


(cosh  rav)b„„  +    (1-km    )a 
^  '    nm        '  '    nm 


^nm    ' 


^^^^^  ^nm  ^"^   ^nm  ^^®  polynomials    In  aj^^^  >\;)  >  >   '^kS    *    \X    ' 
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with  k  <  n,   and    of   weight   n   withrespect  to   the   first    index, 

and   furthermore,    f       =  g        =0     vhenever: 

'  nm  ~  ^nm  - 

m  >  n  for  any  n, 

or     n+m  =  odd  Integer,  any  n  and  m. 

Prom  (14-.7)  ve  obtain 

«       P  2  2 

(Ij.,8)    a   -  a^a   =  f  -(mthmv)g ,  a„  =  m(l-km  )th  m  v  , 

^   '    nm    m  nm    nm  ^     '^nm  '   ra    ^      '       * 

and 

(l;,9)        -   -  ^^"^^ 


J^™   m  sinh  m  v 

Since  the  equations  (l4»8)  are  homogeneous  for  values  of  m 
and  n  listed  above,  and  slnoe  the  initial  data  for  the  a's 
are  also  homogeneous  for  these  values,  one  gets  from  (i+.8) 
and  (U.9), 

(k.lO)   a   =  b   3  0  for  m  >  n,  or  vhen  n+m  =  odd  integer, 
nm    nm  ~  '  ° 

The  third  part  of  the  theorem  implies  then  that  the  coef- 
ficients of  e^  In  (l+.l)  and  (U.2)  are  terminating  series  with 
at  most  [— ]+l  non-indentlcally  vanishing  terms. 

Part  1  of  the  theorem  Is  proven  by  demonstrating  that 
both  Z  and  ^   must  be  even  functions  of  fii    '    Part  2  is  shown 
by  observing  that  incompresalblllty  and  periodicity  o f  the 


(2)  It  is  not  true  however  that  Z(x,t)  vould  be  an  odd  function 
of  X  if  Z(x,o)  were  odd.  Consider,  for  example,  the  same 
problem  and  set  x=(£— ).   Using  E,   as  an  independent  variable 
we  would  be  faced  with  exactly  the  same  problem  with  E,   re- 
placing X  throughout  and  with  the  initial  data  Z(f,o)  = 
=  7(x,t)  =»  e  sin  E   replacing  (2,8«),   The  solution  for  Z  is 

Z(^,t)  -  Z(x,t)  =Z:  e""  H  \   a^^(t)  I(1%1'"^)  cos  m^ 


which  is  odd  in  ^. 
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problem  imply  the   conservation   of   volume    In  a   wave   length. 

Hence 

rZ-n.  CO 

V(t)  =  \   (v-?(x,t))dx  =  constant  =  2it(v~  ^  a  e    )=2nv . 
Jo  n=o 

Consequently,  a   =0  for  all  n.   Proof  of  part  3  preceeds 
by  induction  on  n  and  is  rather  involved.   Since  it  is  simi- 
lar to  that  given  in  an  analogous  case  In  [3]»  Appendix,  we 
omit  it  here. 

In  view  of  the  theorem  only  the  following  coefficients 
are  not  identically  zero: 

^1  '  ^11 

^22  '  ^20  *  ^22 

^31  *  ^33  *  ^31  '  ^33  ®^°* 
As  one  goes  further,  the  n^omber  of  coefficients  to  be  com- 
puted increses  and  also  the  equations  that  determine  them 
become  more  complicated,   i^'hlle  the  work  involved  is  straight 
forward,  the  amount  of  work  in,  say,  fourth  approximation 
is  already  prohibitive.   It  is,  however,  possible  to  char- 
acterise the  results  further  without  an  explicit  determina- 
tion, with  which  we  now  proceed. 

It  turns  out  that  f, ,  =0,  g,,  =0,   Consequently, 
from  equation  (I;, 8),  a, ,  =  cosh  a,  t.   We  use  here  ^Iio 


notation 


a  =  -»-  .pT  ,  a   =  -  fa^  ,  if  a^  >  0 
m    \j'  m '  '   -m    y  m  *     m 


a  =  +i  JTc?T  ,  cc   =  -i  JU^\    ,  if  a^  <  0, 
m      Y'm''-m      \l  '    m'    '  m 
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It  now  can   be   3hot<<n  by   induction  on  n  that   a        Is   made   of 


terms   of  the    form 


P^(t)e^^   . 


Here : 


P      -   Is    a  polynomial    of  order  not  exceeding  n 


a  =  J~   a        ,      with  V"    |m. |    <   n   , 
t"  ^^1  ^        1     - 

In  particular,  .when  all  the  m^^  >  0,  and  ^I  "ij^  =  J^»  one 

shovs  that  P  (t)  la  a  constant.  This  information  now 
n 

suffices  to  determine  the  behavior  of  a   for  large  t.  We 

nm        ^ 

distinguish  two  cases: 

2 
Case  1:  k  <  1,  i.e,  a^  >  0 

We  have , 

m^a?  -  a^  =  m  rm(l-k)thv-(l-m^k) th  m  v]  >mth  m  v  (ni^-.l)k  >  0, 

since  mth  v  >  th  m  v.   Hence 

I  a  I  <  ma, 
'  ra'  —   1 

and  consequently 

1^1  =  IZ  \\  <  a^  21  l"^il  :^  ^^1 

Clearly  then  the  dominant  term  (for  t  large)  in  the  expression 

for  a   is  of  the  form 

na  t 
constant  x  e   1 

This  certainly  implies  divergence  of  the  series  for  t  —>  co  , 

showing  that  either  our  procedure  falls  for  large  t,  or  that 

the  assumption  of  a  well  defined  lower  surface  of  the  layer 

is  untenable  for  large  t.   There  is  no  indication  that  the 

series  (l+.l)  and  (l+.S)  converge  for  any  t.  The  explicit 

results  of  the  next  section  actually  indicate  that  at  best 
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an  asymptotic  behavior  of  these  series  may  be  expected.   Thus 
one  would  not  necessarily  improve  the  accuracy  of  results  by 
Increasing  the  number  of  terms  under  consideration.   While 

we  lack  a  criterion  on  hov  many  terms  It  is  best  to  keep  for 

•J 
each  t,  retention  of  terms  up  to  order  e   leads  to  a  descrip- 
tion in  qualitative  agreement  with  the  observed  behavior  of 

the  layer. 

2 
Case  2:  k  >  1  ,  i.e.  a,  <  0 

2 

Since  in  this  case  a  <  0  for  all  m,  a  is  pure  imaginary  and 

a^   grows  at  most  as  t^»   This  case  leads  to  results  similar 
nm  ^ 

to  those  that  would  be  obtained  if  the  layer  were  accelerated 
In  the  opposite  direction^'^  .   The  latter  leads  to  a  stable 
configuration,  namely,  to  standing  waves.   We  may  therefore 
anticipate  a  similar  situation  in  the  case  k  >  1  ,  that  is, 
when  T  is  large,  or  a  and  I   are  small. 

Further  discussion  will  be  restricted  to  the  case  1 
only. 
5»  Explicit  results  up  to  terms  in  e  « 

Carrying  out  the  program  described  in  "the  preceeding 
section  up  to  terms  in  e   ,  one  gets  the  following: 
Prom  substitution  of  (I;. 1,2)  in  (2.5"'")» 


(3)  In  that  case  the  formulation  of  the  problem  would  be 
essentially  the  same,  except  for  replacing  a  by  |a|  in  for- 
mula for  't   and  k,  and  changing  the  sign  in  front  of  z  in 
equation  (2.6').   On  using  the  same  formalism,  we  would  get 

for  the  a  the  formula  ^         p 

a     =  -m(l+km    )th  m  v    • 
m  ' 

Thus   all   a     VTOuld   be  pure   Imaginary. 
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^22  ~   ^11^11  °°^^  '^ 


f--  =  3a,, boo  cosh  2v+  ^a2pb,,cosh  v-  g-  a,-,b,,    slnh   v. 
Prom  substitution  f\f   (l4..1,2)    in    (2,6' )> 


^gii  =  n 


(5.2) 


§20  =  I  ^iihl^^^   ^-  5  ^11   °°"^^^   -  ^  ^11 


^22  =  I  ^11^1   ^^^^   ^   ■"  ^  ^11 


3   „2 


)^31  =        8  ^^il'"^11^22  ^^^^  2v  -  ^  a^^b^^   oosh  V 

1    •  12 

+  i-ap2b,-]^sinh  '^-t>,,b22  cosh  3v+  -  a^^^b^^  slnh  2  v 


g 


33  =  8  ^^ll''^11^22  ^^^^  2v-  ia^^b^^  cosh  v 


1    • 

+  —  ap2t'-,-,  sinh  v  +  b,,b22  cosh  v  , 

Using  (5.1,2)  In  (U«8)  we  now  get  explicit  equations 

for  the  a   which  can  be  solved  in  terms  of  elementary  tran- 

rim 

cendentals,     VJe  get,   using   the   Initial    conditions    {l].,3,k)t 
a,  ^    =   cosh   a,  t 


1-k 


22 


2  a' 2 
1-k 


I   (cosh  a2t  -   1) 


2  2      2    "1 

(1+ )([|.a,-a2)      (cosh  a2t   -   cosh  2a,  t) 


2  cosh  2v 

(5«3)<(    Sot    =  A   cosh    (a^+a2)t+B  cosh   (ap-a,  )t  +  C   cosh  3ciit 

+  Dt   sinh   a,t   -    (A+B+C)    cosh   a,t 

a-^  =  A   cosh   (a, +ap) t+B   cosh   (cp-a,  )t  +  C   cosh  3ci,t 

+  D   cos    a,t   -    (A+B+C+D)    cosh   c-t. 
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A,B,C ,D,A,B,C ,  and  D  are  numerical  coefficients  (depending 
on  the  parameters  v  and  k)  which  are  too  cumbersome  to  be 
reproduced  here.   Their  exact  values,  however,  were  used  In 
the  calculations  discussed  In  Section  6« 

Simplified  approximate  expressions  for  the  a   can  be 
obtained  for  small  t,  for  large  t,  and  for  particular  values 
of  V  and  k.   For  t  small,  we  get,  using  power  series  expan- 
sions In  t, 


y 


'ii~ 


2 


i^.hK 


^ 


^22^- 


^-^     t^  .ii=k)Ith  V  (U  —2 

2  cosh  2v        12  cosh  2v 

, 2(1-Uk) 


^31- 


33 


(l-k)(l+th^v)cosh  2v 

,2 


if^  th  v(th  V  th  2v  -  1)+  2-  k  th  v 
U.  2   15 

^  {l-k)(3  th  v+th  3v-h|.  th  v  th  2v  th  3v)+  l^-kth  3v  t^ 


16 


16 


For  large  t,  we  get,  keeping  the  dominant  terms 


1  ^1^ 
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(5.5) 


'33 


■»•  th  V  -  th  3v  +  liL-2l__]  +  ap  [cth  2v  ~  th  3v]) 
-  -^  k  th  3v  -  -^afc^-iith  2v  th  3v  +  th  3v  cth  v 


+  2 (3-th  V  th  2v)  ) 


2 

alnh  V 


1  3a  t 


L 


Our  formulae  simplify  considerably  in  the  case  considered  by 
Taylor  [?],  that  is  in  the  case  of  a  sheet  of  infinite  width 
(v  =  0),  neglecting  the  effects  of  surface  tension  (t  =  0, 
hence  k  =  0),   Then  a     =  /  "i"  ,  and 


a,,  =  cosh  t 


/r^    L\     =^n-D "^  X  cosh  ,'?t  -  y  cosh  2  t 

(5.6)^  22     ^   2      >/  U. 

Si       =  -i  cosh^^t  -  ir-  cosh  3t  -  i  t  sinh  t  -  ^ 
31   [i     *     16     ^    8  16 


"5  -5  n  o 

a-_   ~   ""  r  cosh   t   cosh^/^t   +  —  cosh   3t   +  ^^  cosh  t  +  -^  cosh   ^(^"ti 


To   the   order   considered,   the   equation   for   the   profile 
of   the   boundary    is 

2  3 

z  =   Z(x,t)   =    ea, -,   cos    x  +    e   a^p  cos   2x  +   e^(a-,cos    x  +   a,-.cos3x) 

—  K   (cos   X   -   a   cos    2x  +   b   cos    3x). 


Here, 


(5.8), 
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K  =  K(t)   =    tn^^   +   e   a 
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The  shape  of  the  profile  at  any  time  depends  only  on  the 
values  of  a  and  b  and  can  be  followed  by  studying  in  the 
(a,b)  plane  the  locus  (a(t),  b(t)).   Depending  en  the  region 
where  the  representative  point  is  located,  various  shapes 
are  assumed.  We  distinguish  h,   essentially  different  groups 
of  shapes,  described  in t he  table  below  for  half  of  the  wave- 
length, ioe,  for  o  <  X  <  7t,  by  listing  their  consecutive 
minima  and  maxima.   (At  x  =  o  or  at  x  =  ti,  we  have  always 
either  a  minimum  or  a  maximum.)   Since  all  calculations  lead 
to  positive  values  of  a,  we  restrict  our  considerations 
to  a  >  o. 


Table 


Location  of  (a,b) 


Region  I:  lj.a^+(6b-l)^  <  1 
and  -li.a+9b+l  >  0,b  <  2/l5 


Region  II:  I|.a^+(6b-l)^  >  1, 
-[(.a+9b+l  >  0,b  >  ?/l5 


Shape 


max,  min. 


max,  min,  max,  min. 


Region  III:  "hfi,   <  9b+l  <  I|.a    min,  max,  min. 


Region  IV:  l4.a+9b+l  <  0 


i   min,   max,   min,   max. 


See    also   Figure    3» 

In  case  of  sheets    of   finite    thickness    ( v  <   oo  )    a  break- 
down will  occur   at    the   time   when  Z(x,t)    first   reaches    the 
value   V   for  some  x,    since   at    this    time    the  liquid  region  will 
become   disconnected.      At   that   time    the   sheet   will   break  by 
pinching   into   filaments    which  normally  would  soon  acquire 
a  circular  oross-section  because    of  effects    of  surface   ten- 
sion.   (If  initially  Z  had   the    form  Z(x,y,o)=  a  cos  -r  cos  ^ 
the   sheet   would   break   into   spherical   drops.)   Depending  on  the 
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location  of  {a,b)  at  this  time,  the  sheet  m^  break  into  one, 
two,  or  three  filaments  per  wave-length.   Thus,  if  at  that  time 
(a,b)  Is  in  Region  I,  the  break  Is  into  one  filainent  per  wave- 
length;  if  it  is  in  Region  II,  -  into  either  one,  two,  or 
three  filaments, if  it  is  in  Region  III,  -  into  tvo    filaments; 
If  it  is  in  Region  IV  -  Into  two,  or  three  filaments.   These 
statements  presuppose,  of  course,  that  Inclusion  of  terms  up 
to  E-^   only  is  adequate  for  the  description  of  the  sheet  up  to 
the  time  of  breakdoiAn. 
5 .  Discussion  of  results 

A  complete  discussion  of  results  is  very  difficult, 
since  the  problem  contains  three  independent  parameters.   We 
therefore  restrict  our  considerations  to  special,  but  sig- 
nificant cases.    Firstly  we  consider  the  Taylor  case,  v  =  oo  , 
Here,  we  first  choose  k  =  0,  e  =  ,01  in  order  to  compare  our 
results  with  Pennington's  machine  calculations  [6]  who  used 
these  data.   Next  we  chose  e  =  .5  and  a  selection  of  values 
of  k.   To  exhibit  the  behavior  of  not  _too  thin  sheets,  we 
consider  the  cases  v  =  1,  e  =  ,5  and  v  =  ,2,  e  =  ,1,   again 
for  a  selection  of  k.   Lastly  we  consider  the  behavior  of 
very  thin  sheets,  v  «  1  excited  spontaneously,  e  «  v. 
1 .  Taylor  case:  v=OD,k=0,  e=.01 

The  (a,b)  loous  is  shown  in  -figure  i|,  while  the  profile  shapes 
are  drar-n  at  various  times  in  figure  5*   The  dotted  line  is  the 
profile  predicted  by  the  linearized  theory  at  t  =  5.1 
(Z  =  .75  cos  x)  and  shows  that  the  non-linear  theory  predicts 
a  slowdown  of  instability.   Our  results  compare  well  with 
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.2U. 
those  of  Pennington.   In  particular,  the  appearance  of  a 
bulge  near  the  trough  at  t  =  U»8  confirms  the  same  feature 
deduced  from  inactive  computations. 

2,  Taylor  case,  v  =  oo  ,  e  =»  ,S,    k  =  0.  •^.  y^,  l/3,  2/3,  1 

(a,b)  loci  are  in  figure  k,   profile  shapes  in  figure  6, 

3.  V  =  1,  e  =  .5,  k  =  0,  27»  j£»^  2/3,  1 

According  to   the    instability   theory,    this    sheet  will   break 
Into   one  piece    (per   wavelength)    at   t  =   3«0L|.  for  k  =    .1|3,    into 
tvo  pieces   at   t  =   1.1L|.  for  k  =    .1,    and  into   three   pieces    at 
t  =    .66   for  k  =    .037.      The  results    of   our   computations    are 
shoT^n   in  figures    7   and  8,      It    is    seen  that   the   sheet  breaks 
into   one   piece   for  k  =   2/3,    1,    into   two  pieces   for   k  =   0, 
l/l2,    1/3,    ■^nd  into   three  pieces    for  k   =  -^  =    .037.    For 
k  =    .037   the   breakdown   is   at   t  =   1.96,  or  almost   3   times    the 
time  predicted  by   the   instability   theory.      Except   for  the 
case  k  =   0,    these  results   are   thus   roughly   in  agreement   with 
those  predicted  by  the    instability   theory,   except   for   time   of 

breakdown. 

11  12 

1;.   V    =f  .2,    e  =    .1,   k  =   0,  ^,  -j^    »   'j*'^'   ^* 

According    to  the   instability  theory,    this   sheet    vi  11  break 
into  one   piece   at   t   =   13.5   for   k  =    .50,    into   two   pieces   at 
t   =  3*35   for  k=  .13,    into  three   pieces    at   t  =   1.26   for  k  =    .05. 
The  results   of    our  computations    are   shown   In   figures   9   and  10. 
It   is   seen  that   the   sheet  breaks    Into   one  piece   for  k  =  2/3, i, 
tTAO   pieces    for  k  =  •«,    and  into   three  pieces   for  k  =  T^j-^y*   ^* 
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Again  these  results  are  roughly  In  agreement  with  those 

predicted  by  the  instability  theory,  except  for  the  time 

of  breakdovn. 

5.  Spontaneous  disturbances  of  thin  sheets. 

Assume  that  e  is  so  small  that  the  shape  coefficients  remain 

small  until  t  becomes  so  large  that  the  dominant  terms  in 

these  coefficients,  formulae  {5«5),  can  be  used  in  place  of 

the  exact  formulae  (5»3)»   Assume  furthermore  that  v  la  so 

small  that  a  pover  of  v  can  be  neglected  when  compared  to 

a  lover  power  of  v.   It  can  be  shown  that  formulae  (5»U) 

simplify  then  to: 


'11 


lA^ 

■2^ 


1-k   2a,t 
^22  -^  -  l5k^  ^ 


(6.1) 
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^31'"  ' 


l+3k 
5l2kv' 


3a-^t 


I. 


a  r^   (l-k)(5-13k)  J^l^ 
^^    102U  k^v^ 


e 


provided  that  k  =  0  and  k  ^  1« 

For  a,b,  and  K,  formula  (5«8),  we  now  get 

K  =  vk(l-k)c^/a 


(6.2) 


y  a  =  — -i- — 
2-pV 


b  =  -12 

^     2-p2-^ 
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vhere 


(6.3) 


a,  t 
r  ee 


\ 


t 

I 


a  =  1-k 

P^  =  -i  k(l+3k) 


T^ 


(l-k)(5-13k). 


Accordingly,  Z/v  depends  only  on  two  parameters,  namely 
k  and  or  ,  and  this  simplifies  the  discussion. 

For  a  fixed  k,  the  breakdown  of  the  sheet  occurs  for 
that  value  of  a  for  which  max  (Z/v)=l,   The  maximum  of  Z 
occurs  either  at  x  =  0  and  then 
(6.1;) 
or  at 


Z„,o^  =  K(l-a+b) 
max 


(6.5) 

and  then 
(6.6)   Z 


X  =  cos 


-1  a-\/a^-3b+9b^ 


6b 


max 


=  K(- 


27b' 


|b(l+3b)H-lbl 


,_a_.2   1   _1_ 
.6b  ^  h"   12b 


) 


Using  equations  (5)  and  (6),  one  gets  the  following  equations 
for  o"  at  the  time  of  breakdown: 
either 

(6.7)  kr2a-o^  (l-k)+  ^^  o  ^  (5-20k+7k^)  =  1  , 
or 

(6.8)  -8kY(|-)^+k(^)(2+(3Y-P^)o2)+8k|Yl 


(f)2 
L  by 


1_ 
6y 


*<\  *  h-y 


3 
2 


=  1. 


whichever  yields  a  lower  value  of  o"  for  a  given  k.   It  turns 
out  that  for  k  >  .737  equation  (6,7)  must  be  used,  while  for 
k  =  .737  equation  (6.8)  must  be  used.  The  locus  of  (a,b) 
at  the  time  of  breakdown  Is  shown  in  figure  11. 
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The  locus  (a,b)  can  be  obtained  from  the  parametric  repres- 
entation, equation  (6,2),  or  from  an  explicit  equation, 
obtained  by  eliminating  ct,  namely 

(i£  b)2  -  (Me  a)2  =  1 

Y  a 

Per  each  value  of  k  this  locus  consists  of  a  branch  of  a 
hyperbola  passing  through  origin,  A  few  of  these  loci  are 

shown  in  figure  II,   As  k  — ^  0,  the  limiting  position  of 

^  2 
the  locus  is  the  parabola  b  =»  -^a   ,  and  as  k  — ^  1,  the 

o 

locus  shrinks  to  a  point. 

Profiles  at  the  time  of  breakdown  are  shown  in  figure 
12  for  several  values  of  k.   For  k  >  .737  the  sheet  breaks 
into  one  filament  per  wavelength.   For  k  <  .737  we  have 
two  filaments,  except  that  for  very  small  value  of  k(k  <  .5) 
a  breakdown  into  three  pieces  may  occur. 
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